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I. INTRODUCTION 


The development ol a better understanding of the transfer of energy 
from an incoming laser beam to the target material it strikes is & research 
goal of the Wave-Matter Interaction Group at the Polytechnic. To accompany 
our ongoing theoretical research program we have begun an experimental 
program to measure parameters which monitor the state of the target material 
in real time on a nanosecond or faster time scale. Although we would like to 
be able to measure directly parameters such as the target surface tempera¬ 
ture, for example, on such a time scale, it is extremely difficult to measure 
temperature directly even on a microsecond time scale. However, since the 
reflectance of the target can be measured in real time with subnanosecond 
resolution, it may be a crucial and process revealing parameter in following 
the transfer of energy from an intense laser beam to a metal target surface. 

For opaque targets such as metals, the absorbence is one minus the reflectance. 
Thus, the temporal behavior of their surface temperatures, which cannot be 
directly measured on such a fast time scale, can then be calculated from the 
measured reflectance. 

At present, materials processing applications are usually carried out at 

7 2 

incident laser power densities < 10 W/cm to avoid plasma formation at the 
target surface. The effect of a plasma has been unpredictable and difficult 
to control. When located close to the target, the plasma can enhance the 
coupling, and when further removed, it can completely shield the target 
from the laser. With a better understanding of the initiation and energy 
transfer processes during plasma formation, we may expect that materials 
processing applications could utilize the considerably higher power densities, 
as high as 10^W/cm^, which are now available from lasers. Not only could 
applications of increasing scale be considered, but fluid dynamic effects could 
be utilized to make these processes more efficient. 

Our aim then is to relate changes in the target's reflectance to physical 
changes taking place at the metal surface, such as surface deformation, slip, 
melting, plasma formation, ejection of liquid metal, etc. The reflectance 
behavior of a metal surface irradiated by an intense laser pulse must be 
consistent with the laser energy absorbed, the metal surface temperature 
attained, and the surface damage produced. 
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During this AFOSR program, we have measured both the specular and 
total reflectances of vapor'deposited copper surfaces during irradiation by 
a Q-switched ruby laser. We have calculated the las sr energy absorbed and 
determined the temperature of the copper target surface using a one-dimen¬ 
sional heat conduction approach. Finally, we have related the reflectance and 
temperature histories to any permanent change or damage observed in the 
copper target material. 











II. TARGET REFLECTANCE BEHAVIOR 


Experimentally, the reflectance of a metal surface has been reported 

to undergo a sharp and substantial decrease during an intense laser pulse. 

For incident laser power densities above threshold values which lie between 
7 9 2 

10 and 10 W/cm , the measured reflectance of a metal surface drops sig¬ 
nificantly before partially or totally recovering as indicated by the curves in 
Figure la. This type of reflectance change was first repeated by Bonch- 
Bruevich et al. * during individual spikes •> 400-ns FWHM (Full Width at Half- 
Maximum intensity) within a millisecond 1.06-iiNd glass laser pulse on targets 
of silver, copper, aluminum, dural and steel as indicated. 


There are three distinct regions within a generalized reflectance curve 
as illustrated in Figure lb: first, AB - an initial steep decrea.-ie to - one-half 
of the initial reflectance; second, BC or BC 1 - a plateau region during which 
the reflectance remains approximately constant; and third, C'D' - a complete 
reflectance recovery, or for higher intensity laser pulses CD, a further de¬ 
crease to - one-tenth of the initial reflectance followed by a partial recovery. 

4 

Zavecz, Saifi and Notis reported similar reflectance curves during 
Q-switched 1.06-nNd:YAG laser pulses - 60-ns FWHM on single crystal 
copper and tantalum targets as shown in Figure lc. The remarkably close 
agreement between Zavecz's 0, 64 plateau reflectance value znd the 0. 66 



Time (ps) 

(a) Reflectance 
measurements of 
Bonch-Bruevich et al. 



Time 

(b) General reflectance 
curve. 



(c) Reflectance 
measurement of 
Zavecz et at. 


Fig. 1. Metal reflectauce behavior during an intense laser pulse. 
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value of Bonch-Bruevich for copper, the only metal common to both investi¬ 
gations, is evident in Figure 1. Plateau features have also been reported by 
a 

Dymshits, using a 30 -ns FWHM 1.06-|i Ntl glass single laser pulse on thin 
aluminum films, and Walters and Clauer, using gain-switched 50 to 70 - ns 
10. 6-p TEA COj laser pulses on polished poly crystalline aluminum rods. 

Bonch-Bruevich et al. * associated the initial steep decrease in the re¬ 
flectance curve, AB, with the heating of the metal surface to melting and 
with the drop in electrical conductivity which accompanies melting. The 
plateau BC or BC', they suggested, indicated that the temperature of the 
molten layer remained constant while the absorbed radiant energy propagated 
a melting wave into the solid. The further reflectance drop during region 
CO was attributed to a decrease in the amount of energy conducted to the 
liquefaction wavefront and, therefore, indicated a second region of increasing 
surface temperature. 

4 

Zavecz et al. tested Bonch-Bruevich's conjecture that the metal was 
melting during region BC. A single-crystal copper target was irradiated 
with a laser intensity (3.2x 10 W/cm ) which was sufficient to attain the 
plateau reflectance BC, as shown in Figure lc, but not high enough to cause 
the second reflectance decrease. Examination of the target by scanning 
electron beam microscopy revealed no evidence of melting. Melting occurred 
during the second decrease in reflectance, Ct of Figure lb. A major dif¬ 
ference between the two experiments was the use of an integrating sphere by 
Bonch-Bruevich and, therefore, a total reflectance measurement while Zavecz 
et al. measured only the specular reflectance. In ensuing experiments to 
Zavecz's, Koo and Slusher^ proposed that the reflected beam could have been 
deflected into a wide cone of angles by a transient grating structure produced 
on the metal target surface. Ready, ^ using a 10.6-ji CO, laser, and von 

Q “ 

Allmen using a 1.06-p. Nd laser on metal targets, found substantial and 
permanent changes from specular to diffuse reflectance during the laser pulse, 
Ready used a probe laser to subtract plasma absorption and concluded that 
no recovery of reflectance occurred during or after the pulse. An ellipsoidal 
light collector was employed by von Allmen in front of a copper target in air. 
Ke found no significant change in the total reflectance until the metal surface 
exceeded the boiling point. Prokhorov et al. ^ suggested the formation of a 
liquid dielectric layer causing an enhanced absorption. More recently, 








8 

Dymshits observed Figure I-type reflectance changes during a specular re¬ 
flectance measurement with the metal target film in vacuum. He suggested 
that the reflection takes place from a plasma formation in front of the metal 
surface. Walters and Clauer 10 report a sharp decrease to -0.35 in specular 
reflectance at the melting point of AI followed by nearly complete recovery. 
No luminosity was observed and no significant permanent loss of the high 
specular reflectivity was detected. 

An under standing of the metal surface reflectance behavior during an 
intense lasei pulse is lacking at present. In addition to Bonch-Bruevich's 
suggestion that the initial steep decrease in reflectance is associated with 
the heating of the metal surface to melting, three classes of alternate ex¬ 
planations have been offered for the steep decrease in metal target reflec¬ 
tance during an intense laser pulse: 1) surface deformation, 2 ) plasma 
formation, or 3) nonlinear process. Before discussing these, we shall 
first examine the change in reflectance expected when a metal is heated to 
its melting point. 






1H. FREE ELECTRON MODEL 


The Fresnel expression for the reflectivity, R, of radiant intensity at 

12 13 

normal incidence on the surface of a homogeneous substance is * 


r . iiii 

( 5 + 1 | 


2 
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where n is the complex refractive index. From the viewpoint of optics, the 
complex index of re fraction may be expressed in terms of the optical constants 
n and k of the material 

n = n - ik 

where n is the index of refraction and k is the extinction coefficient. 

In the Drude or free electron model, a metal is viewed as a gas of free 
electrons which interacts with a background positive-ion lattice represented 
by phonons. The complex index of refraction then may be expressed in terms 
of the plasma frequency 

tu p * J 4rrNe^/m* 


the average electron-phonon collision frequency a \> c and the frequency of 
the incident light (b,^^ s 2nc/X, 


tb / v_ \ 

"- 1 -rr— E —r(‘ + i i—:•)* 1 - 


"2 - 1 

"light + v c 


light 


1 ,, . . v, 

rT (1 + l i ) 


<b” + v 


where tb s tu ../<M , and v a V /tb . Here N is the free electron density, 

light' P * c P 

e is the electronic charge, m is \he effective mass of the free electron in 
the metallic lattice and i a . For visible light incident on a good con¬ 
ductor 


tb * 10 16 sec” 1 > tb. * 10 15 sec" 1 > v «10 14 sec -1 . 
p laser c 

Therefore, tb - 10’* and V - 10’^. 

When the temperature is raised to the melting point, no significant 
change in the free electron density is expected for a metal whose ban&gap 
is much greater than the thermal energy change. Therefore, no appreciable 






change is expected in the plasma frequency The phonon density, however, 
may increase significantly and hence also may the electron-phonon collision 
frequency and v, Figure 2 shows the dependence of the Fresnel reflectivity 
on the collision frequency for neodymium and ruby laser frequencies normally 
incident on a copper surface. 

There has been considerable variation in the room-temperature re¬ 
flectance measurements of many metals. Since the absorption depth, X/2rrk, 
is only a few hundred angstroms at most for metals in the visible-near infra¬ 
red spectral region, the measured reflectances of metals are very sensitive 
to the method of preparation of the surface as well as to the history of the 
surface from preparation to measurement. When the metallic surface is 

produced by polishing, mechanical stress can deform the surface by producing 

14 

an amorphous layer which may be several hundred angstroms thick. Even 
more detrimental is the embedding of polishing-compound particles which is 
more likely to take place in soft metals: copper, aluminum, etc. If an oxide 
layer is allowed to form, the reflectance measured will not be that of the pure 
metal. For metals which are good conductors and therefore also possess 
high reflectivities, processes which modify the metal-air or metal-vacuum 
interface generally reduce the reflectance. The highest reflectances of 
metals have been measured for evaporated coatings prepared in a high vac¬ 
uum at a high rate of deposition of the metal and measured immediately after 

15 

preparation while 3till in a high vacuum. Reflectances of Ag, Cu and Af 
measured under these conditions at 1.06|i are .994, .985, and .944, re¬ 
spectively. 

Figure 2 indicates that for a room-temperature copper reflectivity of 
0,985 at 1.06n, the ratio of the electron-phonon collision frequency to the 
plasma frequency v= 0.00749. According to the free electron model then, 
for the reflectivity of copper to decrease to 75%, for example, v must in¬ 
crease so that v> 0.15, This means an increase of more than 20 times in 
the electron-phonon collision frequency. 

13 

Ujihara has examined the temperature dependence of a Drude-type 
free-electron model for high conductivity metals. He asked whether the in¬ 
crease with temperature of the electron-phonon collision frequency can ac¬ 
count for a reflectance decrease to the plateau value by the time that the 
metal reaches its melting point. Ujihara assumed a Debye model for the 







Fig. 2. Reflectivity of copper as a function of the 
electron-phonon collision frequency v c for 
two incident laser frequencies u>. 


phonon spectrum and N-scattering processes on a spherical Fermi surface. 
Both the Debye temperature and the electron distribution were considered 
to be independent of temperature in the 300°K to melting point region. 
Ujihara's results for several metals at ruby and Nd wavelengths are shown 
in Figure 3. The reflectance of solid copper at 1.06|i decreases from 0,95 
at room temperature to 0.73 at the melting point, 1356°K. Since the con¬ 
ductivity of copper, like many other metals, decreases by a factor of 2 upon 
melting, ^ a reflectance of 0. 58 would be expected for liquid copper at 
1.06(1 at the melting point. These reflectance values bracket the plateau 
reflectances of 0. 66 and 0. 64 observed by Bonch-Bruevich and Zavecz, and 
suggest that the Drude-type free-electron model may be able to explain the 

7 

reflectance behavior observed. Indeed, as shown in Figure 4, Chan et al. 
claim to have confirmed Ujihara’s theoretical reflectivity curves using a 
ruby laser on copper and aluminum targets and calculating the surface 
temperature on the basis of a one-dimensional heat flow model. 
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Fig. 3. Ujihara's calculated re¬ 
flectivities (13) for 
several metals at two 
laser frequencies. 


Fig. 4. Chan's measured reflectances 
(7) of copper (triangles) and 
aluminum (dots) at 0.69|i. 


17 

Ujihara's analysis has been examined by Walter, who pointed out 
that Ujihara, in his determination of the plasma frequency U) , used the 

18 P 

1913 effective-mass measurements. For copper, the old value is almost 

two times the more recent measurements. When the resulting low value for 
the plasma frequency was used with values of the optical constants, n and k, 
the room-temperature value Ujihara calculated for the electron-phonon col¬ 
lision frequency was 2. 5 times too high and the room temperature reflectivity 
at 1.06|jl was 95.2% instead of 98.5%. With these incorrect initial values, 
his analysis then predicted substantial changes in the collision frequency and 
in the reflectivity which are not confirmed by experiment, as indicated in 
Figure 5. 

The solid curve in Figure 5 shows the temperature dependence of the 

reflectivity of copper calculated by Ujihara. The dotted curve extends 

Ujihara's analysis down to 77°K, The optical properties of copper have 

_g 

been measured polarime* rxcally in ultra-high vacuum (-10 Torr) by Pells 
and Shiga^ over the temperature range 77 to 920°K. Reflectances com¬ 
puted from Pells and Shiga's experimental values of n and k are also indicated 
in Figure 5. The discrepancy between theory and experiment is substantial. 
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However, when the improved data listed in Table 1 is used, then the theory 

of Ujihara is in much better agreement with Pells and Shiga's experimental 

results. The dominant change responsible for most of the improvement is 

20 

Schulz's careful determination of the effective electron mass in copper 


Table 1. Input 

Data for Ujihara's 

Theory 

Parameter 


Original 

Ujihara 

Analysis 

Improved 

Data of 
Reference 17 

Index of refraction 

n 

0.150 a 

0. 1074 b 

Extinction coefficient 

k 

4.049 a 

3.9104 b 

Wavelength 

X 

0.70 |i a 

0. 65 |i b 

Effective mass 

* 

B 

2.56 C 

1.45 d 



a Reference 21, 
“Reference 22. 
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c Reference 18 
^Reference 20 











including a correction for the anomalous skin effect. With these improve¬ 
ments in the input data used in the Ujihara theory, we may conclude: 

23 

(1) the reflectance of liquid copper is -85%, not -65%, and 

(2) it is unlikely that the temperature dependence of a Drude-type 
free-electron model can explain the substantial reflectance 

17 

changes reported for high-conductivity metals such as copper 


In Figures 6a and 6b, the temperature dependence of the optical con¬ 
stants n and k, and the free-electron model constants Ui and v , are dis- 

P c 

played and compared with the polarimetric measurements of Pells and 
19 

Shiga. As expected, it is evident that the plasma frequency does not 
change significantly with temperature. The plasma frequency can be de¬ 
termined from the optical constants n and k or from the free electron density 

* 22 
and effective mass N and m . When the most precise n and k values are 


used to determine U) for Ag, Cu and Au, good agreement is obtained with 

P 

U)p values determined from N and m , using the effective mass values of 
Schulz. Also, the effective mass for electrons in Cu determined by Schulz 


from optical properties at 2(.i is nearly identical to the value from electronic 

20 

specific heat determinations. It is also evident in Figure 6b that the 
electron-phonon collision frequency undergoes a much smaller increase (-4) 
when copper is heated from room temperatxire to its 1356°K melting point 


than that required (> 20) to produce the substantial reflectance drop to 


Figure 1-type plateau values. 


Finally, it is clear from the plots of n and \> c in Figure 6 that even 
with the improved input data listed in Table 1, the theory of Ujihara does 
not have the correct functional dependence with temperature for either the 
index of refraction or the electron-phonon collision frequency for T< 500°K. 
Refinements such as a departure from sphericity of the Fermi surface and 
more than one electron-phonon collision frequency could be considered to 
improve the Ujihara treatment. 
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Fig. 6. Temperature dependence in copper cf (a) optical material proper¬ 
ties: index of refraction n and extinction coefficient k, and (b) free- 
electron material properties: plasma frequency u and electron- 
phonon collision frequency v . P 
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IV. REFLECTANCE MEASUREMENTS 


Other explanations for the reflectance behavior which is shown in 
Figure 1 fall into three categories: (1) a surface deformation which direct? 
the reflected beam away from the detector, (2) plasma formation which ab¬ 
sorbs or scatters the reflected light, or (3) a nonlinear process causing en¬ 
hanced absorption within the metal. The first two categories involve pro¬ 
cesses which produce a reduction in light reaching the detector that could be 
misinterpreted as a decrease in the reflectance of the solid target surface. 
The importance of the first category, surface deformation, may be evaluated 
by comparing specular reflectance with total reflectance. The initial ex¬ 
periments of Bonch-Bruevich (Figure la) apparently measured total re¬ 
flectance while those of Zavecz (Figure lc) were specular reflectance mea¬ 
surements. 

Experimentally, the reflectance of the metal surface is obtained by 
dividing the reflected pulse by the incident laser pulse. For a 30-ns laser 
pulse, even a time registration error as small as 1 ns between these two 
pulses can cause a significant error in che measured reflectance. This is 
indicated in Figure 7 where the solid pulse curve represents 30-ns FWHM 
Gaussian incident and reflected pulses, and the solid horizontal line is the 
corresponding reflectance. Time registration errors of 1 and 3 ns are 
shown by the dashed curves with the dashed reflected pulses now preceding 
the solid incident laser pulse. A shift of only l ns, corresponding to 3% of 
the width, can cause a reflectance error as high as 25% in the wings of the 
pulse shape. 

To avoid a time registration error, we developed a mc-thod which uses 
a single detector- for both pulses. A sample of the incident laser pulse is 
sent through an air delay path. ** Both specular and total reflectance mea¬ 
surements have been carried out on copper targets. The experimental ar¬ 
rangement is sketched in Figure 8. A TRG model 104 ruby laser operating 
in the Q-switched mode produced 0. 6J in a 30-ns FWHM pulse. A one to 
four beam expander reduced the measured beam divergence to -1 mrad. 

A high-intensity beamsplitter (Melles Griot Cat. No. 03 BTF 007 with 

o 

60% transmission at 6943 A) provided a sample of the incident pulse. This 
sample pulse was delayed through an air path of - 20m for specular re¬ 
flectance measurements and ~ 50m for total reflectance measurements. 
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reflected and incident laser pulses with 0, 1, 
and 3-ns time registration errors. 



(a) For specular reflectance 
measurement. 


Air path 
delay 



attenuator Vacuum chamber 
W * window 


(b) For total reflectance 
measurement. 


Fig. 8. Experimental arrangement. M = mirror, L = lens, F = filter or 
attenuator and W = window. 


- 14 - 











Both the reflected pulse from the target and the delayed sample incident 
laser pulse were detected by the same detector, a TRB 105B vacuum-photodiode 
with a detector area of 5 cm diam, a risetime < 0.3 ns and a dc-1.25 GHz 
bandwidth. The detector output was coupled to a Tektronix 519 oscilloscope 
which has a dc-1 GHz bandwidth. Therefore, the resolution of the detected 
light signal was -0.5ns. 

The copper target surfaces were produced by vacuum-deposition on a 

glass plate of good optical quality. Evaporation and rf sputtering have each 
, o 

been used to produce copper films - 2000 A thick. T'.e sputtered copper 
deposit exhibited stronger adhesion to the glass base* plate than did the 
evaporated coating. Following vacuum deposition, the copper deposit on the 
glass plate was backed up by electroplating to produce a mechanic ally-strong 
copper target a few mm thick. Thus prepared, the copper target was kept 
attached to the glass plate until required for the reflectance measurements. 

Then the glass was detached and the vacuum-deposited copper target placed 
in the vacuum chamber. In this way an oxide-free coppnr surface could be 
prepared and protected until the measurements were taken. As indicated in 
Figure 8, the targets were located in a vacuum chamber, and a mechanical 
pump reduced the chamber pressure below 0. 1 Torr to avoid air breakdown 
at the target surface. 

The specular reflectance measurements are indicated in Figure 8a. 

The 40% reflected beam from the beamsplitter was focused onto the target 

with a 50-cm focal length 15-cm diam quartz lens. At the focused spot -0.5- 

9 2 

mm diam a maximum intensity of 10 W/cm was obtained without any filters. 
This was determined from a measured energy of 75mJ on the target. Corning 
glass filters were used as attenuators to adjust incident reflected and detected 
intensities. Four aluminized mirrors were used to provide an air path delay 
of - 65ns for the sample of the incident laser pulse. 

In the total reflectance measurements, an integrating sphere was used 
as indicated in Figure 8b. The sphere was fabricated from two plastic hemi¬ 
spheres of 20 cm i.d, which were coated with Eastman White Reflectance 
24 

Coating. The Eastman Coating, a composition of barium suphate, binder 

and solvent, is characterized by high reflectance, high stability and nearly- 

25 

perfect diffuse reflection. The manufacturer's specifications show .992 

o 

reflectance at 7000 A, Not only does an integrating sphere carry out a 
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spatial integration of the light reflected from a target surface, but it delays 

the arrival of light at the detector by varying times depending on the path 

within the sphere. Thus the temporal shape of the detected light signal may 

be significantly modified. A 33-ns FWHM ruby laser pulse reflected from a 

metal target surface onto the 20 cm diam integrating sphere was lengthened 

by - 10 ns when measured by the fast vacuum photodiode located at the top of 

the sphere. The pulse lengthening is in good agreement with the theory of 

2 6 

the integrating sphere for pulsed light sources. Both the reflected pulse 
and the sample laser pulse were directed into the integrating sphere for the 
total reflectance measurements. Because of the pulse-lengthening property 
of the integrating sphere, a longer air delay path, - 50m, was rea- red to 
avoid overlapping the reflected and sample laser pulses. After tt t air path 
delay, as shown in Figure 8b, the sample of the incident laser pulse was 
sent through a second entrance window directly onto the integrating sphere 
surface; however, it was not focused on the sphere. 

For each reflectivity measurement, a Polaroid picture was taken of 
the Tektronix 519 oscilloscope trace showing a reflected pulse followed by a 
sample incident laser pulse. To facilitate data analysis the pictures were 
enlarged by a linear factor of - 6, as shown in Figure 9. 


Fig. 9. Enlargement of an oscilloscope trace showing the specularly-reflected 
pulse from an evaporated-copper target surface and a sample of the 
incident 3.15x 10^W/cm^ ruby laser pulse. The horizontal time 
scale is 50 ns/large division. 
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Readings were taken at 2mm intervals of both the upper and lower 
boundaries of the thickened curve trace. Approximately 90 upper and lower 
curve trace values were entered into a PDP 11/60 minicomputer as two 
vectors. An interactive graphic language was then used to average the two 
vectors, remove the time delay of the sample laser pulse, divide the re¬ 
flected pulse by the sample laser pulse, apply the calibrated scale factor to 
change from relative to absolute reflectance and, finally, disnlay the reflec¬ 
tance curve. The calibration factor involved filter, window and lens trans¬ 
missions, mirror reflectivities in the air delay path and the beamsplitter 
transmission and reflection. For incident laser pulses within a narrow range 
of intensity, the reflectance curves were averaged and examples are pre¬ 
sented in Figure 10 for specular reflectance and in Figure 11 for total re¬ 
flectance of a vapor-deposited copper surface for three incident laser power 



Time (ns) Time (ns) Time (ns) 


(a) 2.37x10 a W/cm 2 Cb) 5.56x10 8 W/cm 2 Cc) 9.0x10 8 W/cm 2 

Fig. 10. Specular reflectance of copper. 



Time (ns) Time (ns) Time (ns) 


(a) 2.02x10 8 W/cm 2 (b) 5.64x10®W/cm 2 Cc) 8.49x10 8 W/cm 2 

Fig, 11. Total reflectance of copper. 
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densities, *2, 5, and 9x10 watts/cm . The curves shown are a normalised 
incident laser pulse, the reflected pulse and the metal target's reflectance 
behavior which is the ratio of the two pulse curves. The laser energy ab¬ 
sorbed was calculated from these results and the temperature history of the 
metal surface determined using a one-dimensional heat-conduction approach 
as described in Section V. 







TEMPERATURE OF THE METAL SURFACE 


When laser radiation strikes a metal surface, the light not reflocted is 

absorbed in a thin layer 10”®- 10 cm thick. The energy is absorbed by 

,o 

electrons which equilibrate with the lattice vibrations in a time << 10 sec. 
For nanosecond or longer laser pulses, therefore, the laser radiation ab¬ 
sorbed may be considered to be immediately converted into heat and available 
for conduction inti the metal. For a laser pulse of ~ 30ns duration, the thick¬ 
ness of the region heated during the laser pulse is * jKt or 10"* cm which is 
very small compared with the thickness of the usual metal target, (k is the 
thermal diffusivity, 0.91 cm /sec for copper.) Therefore, the absorbed 
energy may be treated as a surface heat source propagating into a metal of 
semi-infinite dimensions. The heat flow equation is: 


3T(r,t) 2 *q o (£..t) .a .7 

= k v T(r, t) + - 4 - ~ - e- a * 


“St" 


pc. 


( 1 ) 


where T is the temperature of the metal as a function of position on the surface 
£and time t; k is the diffusion coefficient, p the density of the metal, c^ the 
specific heat of rnetal at constant pressure, Qt the attenuation coefficient and 
q o (r,t) is the intensity of the absorbed laser radiation on the surface. The 
dimension z is the moving coordinate into the metal perpendicular to the 
surface, and z=0 is the boundary surface on which laser radiation is incident. 
The boundary conditions are dT/dz = 0 for z < 0, that is, no return heat flow, 
and T(z -* oo,t) = T q . The initial condition is T(r^, 0) = T Q where T q i3 the 
ambient temperature. We shall assume a Gaussian spatial intensity distribu¬ 
tion of the absorbed laser radiation 



e 



2 

» q(t) e 



(?) 


where P(t) is the power absorbed from the incident radiation, w is the radius 
of the irradiated spot and r is the distance from the beam axis on the surface 
plane. Since the attenuation coefficient is very large (Cl » 10^- lO^cm" 1 ), or 
equivalently, the skin depth is very small, the absorbed energy can 

be regarded as a surface heat source, i.e., Ole -» 6(z) for large Cl. In 
cylindrical coordinates (r, <p, z), Eq. (1) becomes 
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with the boundary and initial conditions already given. 

The Green's function for the heat flow equation in a setui-infinite con¬ 
ductor with a surface source having a spatial Gaussian heat distribution, i.e., 
the solution of the heat conduction equation for a unit input of the form 


1 

—4-e w 6(t)5(z) 
TTW 


is shown in References 27 and 28 to be 

1 


z 2 

4<ct 


G(r,*,t) = —372-——*- 

TT pCpA//Ct(w^+ 4/Ct) 


T7* 

w + 4#ct 


(4) 


(5) 


The solution of EJq. (3) for absorbed power P(t) is given by the time convolu¬ 
tion of the Green's function and P(t). 


T(r, z, t) a T q + / G(r, z, t-1*) P(t # ) dt* . 


( 6 ) 


Thus, 


T(r, z, t) s T + -L-/ -i 


z 2 r 2 

.2 " 4*(t-t') ‘ ~Z 




V^rf o pc 


w +4*(t-t')dt 

C» 

(7) 


where 


4(4 = -PCI- 

ttw rr [(w + 4jc(t-t / ) 3 


3ince w» V4ict for w« 0. 1 cm, #c * 1 cm 2 /sec and t - 100nsec. The tem¬ 
perature of the metal surface T(t) = T(0, 0,t) is 


t 

T(t) = T + r 

u ^ 

o 


1 

./rrpc iC 
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dt # 


( 8 ) 


- 20 - 












where K® pc^K is the thermal conductivity of the metal. In the presence of 
energy loss mechanisms, such as melting, vaporization, etc,, q(t') is re¬ 
placed by a net heat input f(t') * q(t *)~i (t*). The loss term i(t) represents 
all heat loss processes besides conduction and is treated as a function of the 
surface temperature, i.e., l[T(t)]. For convenience, a Gaussian spatial 
loss distribution has been assumed. Thus, the general expression for T(t) 
is given by 


T(t) * T + j 


1 - ‘iild.’ 


° o JPV 

Not only is the solution given by Eq. (9) more compact than a LaPlace trans¬ 
form solution, but Eq. (9) also has a clear physical interpretation since 
f(t),/itt represents the energy absorbed in the affected volume and l/pc^ is 
the specific heat factor which converts the absorbed energy to a temperature 
change. Moreover, the integration in Eq. (9) can easily be carried out for 
any function f(t) provided the function f(t) can be approximated by a polynomial 
expansion in t for each region of the integration interval. 

As the incident laser radiation is absorbed, the surface temperature 
of the metal increases until it reaches the melting point. Then while the 
metal absorbs the energy corresponding to the latent heat of fusion the sur¬ 
face temperature remains constant. The thickness of the molten layer de¬ 
pends on the intensity of the absorbed laser radiation and on the thermal 
properties of the metal. After a molten surface layer is formed, the tem¬ 
perature of the liquid metal increases further until substantial vaporization 
occurs. An appreciable amount of the absorbed energy is then spent in 
supplying the evaporation energy and excess kinetic energy of the vapor which 
significantly reduces the surface heating rate. Unless the vapor has a low 
ionization energy, the vapor cloud remains transparent to the incident beam. 
On the other hand, the resultant high vapor pressure exerts a force on the 
molten layer under the irradiated spot. When the tangential shear force due 
to the high vapor pressure exceeds the surface tension of the liquid metal, 
the molten layer starts to move out of the irradiated zone. At this point 
perhaps it should be noted that the effects of radiation pressure of the inci- 
dent beam for the laser pulses considered (10 - 10 W/cm ; »30nsec FWHM), 
and the surface tension gradients on the molten layer arising from the non- 
uniform temperature distribution due to the Gaussian laser-beam intensity 
distribution assumed across the irradiated region are smaller by many 
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orders of magnitude than the effects of high vapor pressure. 

When liquid metal is removed from the irradiated zone, part of the ab- 
sorbed laser energy is required to melt and heat the surface of the freshly 
exposed spot to the temperature of the departed surface layer. From the 
point of view of an observer in a coordinate system located on the metal 
surface, the ejection of liquid metal slows the conductive heating process. 

When the intensity of the absorbed laser radiation is high enough 
(>10^W/cm^) for the metal surface temperature to increase beyond the 
boiling point of the metal, the incoming laser radiation may interact with 
the dense vapor cloud in front of the surface. Now the reflection of the 
incident beam takes place in a region with a gradual change of density, and 
consequently of refractive index, instead of at the sharp boundary between 
two distinct phases. At this point the absorption mechanism of the incident 
laser radiation is switched from the metal surface to the dense vapor layer 
xn front of the metal. A further increase of the vapor temperature can re- 
suit in appreciable excitation, ionization, and plasma formation. The plasma 
can now increase the coupling between the laser radiation and the metal 
surface. Later, when the plasma moves away from the metal target surface, 
a decoupling will take place. 

29 

Temperature histories were determined for the metal targets using 
Eq. (9). The laser power absorbed, q(t), was determined from the measured 
incident laser pulse and the measured total reflectance (e.g. Fig. 11). Heat- 
loss processes in addition to conduction, i (t), were treated as a function of 
temperature and terms were included for melting, vaporization, expulsion of 
liquid metal and solidification in the appropriate temperature regions. The 
thermal conductivity and specific heat were also changed from measured 
values for the solid to those for the liquid in the appropriate temperature 
region. The integration of Eq. (9) was carried out using a trapezoidal ap¬ 
proximation and an iterated substitution method. Table 2 gives the tem¬ 
perature histories for the three ruby laser - copper target interactions whose 
reflectances were measured in Figure 11. Surface temperatures above the 
boiling point of copper, 2840°K, are probably excessive because additional 
loss mechanisms such as absorption by the dense vapor layer and plasma 
formation are not included in this calculation. 
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Table 1. Calculated surface temperature at the center 
of the interaction site irradiated by a ruby 
laser pulse producing 2. 02, 5. 64 and 8. 49 x 
10 8 W /cm'* on a copper target. 


TIME 


SURFACE TEMPERATURE 


’nsec] 

2. 02 MW/cm 

5. 64 MW/cm 2 

1 8.49 MW 

0 

300°K 

300°K 

300°K 

5 

312°K 

343°K 

346°K 

10 

333°K 

483 °K 

464°K 

15 

378°K 

784°K 

910°K 

20 

475°K 

1236°K 

1577°K 

25 

640°K 

2373°K 

4841°K 

30 

811°K 

3730°K 

7869°K 

35 

955°K 

4247°K 

8697°K 

40 

1040°K 

4255°K 

8289°K 

45 

1040°K 

3966°K 

7 484°K 

50 

1012°K 

3606°K 

663 6°K 

55 

971°K 

3368°K 

5849°K 

60 

912°K 

3202°K 

5227°K 

65 

• 

851°K 

3037°K 

4743°K 

• 

• 

130 

572°K 


2352°K 



















VI. DISCUSSION 


The measured total reflectances of evaporated-copper target surfaces 
3how a slow decrease with increasing incident laser intensity or increasing 
surface temperature up to the boiling point of the metal. At low incident 
laser power densities (i.e., < 3x10 W/cm ) the specular and total reflec¬ 
tances of the metal target surface are similar (cf. Figure 10a with Figure 
11a). The one-dimensional heat-conduction calculations indicate that the 

copper surface did not reach the melting point with this incident laser power 

8 2 

density; namely, 2x10 w/cm . No permanent damage was observed in the 
copper. The copper target was examined under 200x magnification and no 
trace of melting was observed. The lack of any permanent observable damage 
is consistent with the calculated maximum surface temperature of 1040°K 
and with the equality of the specular and total reflectances. 

8 2 

At an incident laser power density of 5. 6x 10 W/cm , the specular 
reflectance undergoes a sharp drop which is not observed in the total reflec¬ 
tance (cf. Figure 10b with Figure lib). The one-dimensional heat-conduction 
calculation indicates that the maximum temperature attained by the copper 
surface was in excess of the boiling point. Permanent damage was observed 
at the interaction site. Once again the permanent changes observed in the 
target are consistent with the surface temperatures calculated and the re¬ 
flectances measured which reveal a decrease in the specular accompanied 
by an increase in the diffuse component of the reflectance, 

8 2 

An additional increase in incident laser power density to -9x10 W/cm 
continues the decrease in specular reflectance (Figure 10c compared with 10b), 
Permanent target surface damage is more severe and the crater site is more 
extensive. The recovery in the total reflectance behavior which is evident in 
Figure 11* is caused by a "light flash" produced by emission of light from a 
plasma that is formed at or in front of the metal surface. The spectral con¬ 
tent of this light includes sharp emission lines from neutral copper atoms, 
copper ions and copper dimers (Cu^ molecules). 

When an interference filter is inserted in front of the TRB 105B vacuum 

o o 

photodiode thereby restricting the light detected to i 50 A around the 6943 A 
ruby laser wavelength, then the recovery in total reflectance behavior so 
evident in ure 11c and just noticeable in Figure lib disappears, ^ 
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VII. PROBE LASER REFLECTANCE MEASUREMENTS 


When light from the high-power laser pulse is used to provide the 
optical measurement of target reflectance, the signal is only present and 
hence the reflectance measurement can only be carried out during the short 
time period of the laser pulse length. During the initial portion of the rise¬ 
time and the final portion of the trailing edge of the laser pulse, the signal- 
to-noise ratio of the reflectance measurement deteriorates. Thus, precise 
reflectance measurement can only be made for a time period somewhat 
greater than the FWHM which was - 30ns for the ruby laser used in our ex¬ 
periments. If a cw probe laser were focused to a spot on the target which 
is entirely within the focused spot of the ruby laser, then a second set of 
reflectance measurements could be taken using the cw probe laser. Probe 
laser reflectance measurements would provide two types of additional in¬ 
formation to assist in our understanding of the processes involved: 1) re¬ 
flectance measurements over much longer time periods, and 2) reflectance 
measurements at a second wavelength. Reflectance measurements over 
longer time periods would provide more sensitive measurements during the 
initial stage when the high-power laser is turning on as well as through the 
maximum target surface temperature period as the high-power laser is 
turning off and also into the cooling region as the surface resolidifies. Re¬ 
flectance measurements at a second wavelength can also be important in 
identifying the role of nonlinear processes in accounting for the steep reflec¬ 
tance drop shown in Figure 1, Both reflectances would be expected to be the 
same except when a nonlinear process such as enhanced absorption causes 
a decrease in reflectance at the high-power laser frequency but not necessarily 
also at a second, non-resonant low-power probe laser frequency, 
o 

The 6328 A helium-neon laser is a readily available cw probe laser. It 

o 

is close enough to the 6943 A ruby laser to produce nearly the same response 
from our detection system and is also close enough to avoid exciting inter¬ 
band transitions in targets like copper. In addition, the visible output 

facilitates the alignment of its focused spot to ensure that it is entirely 

o 

within that of the ruby high-power laser. The 615 A separation between 

the two laser wavelengths is sufficient for easy discrimination by ordinary 
o 

100 A commercial interference filters. To separate the signal from a 1 mW 
He-Ne laser from a 10MW ruby laser pulse, we had expected to use a stack 
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of four interference filters since the skirt-to-peak transmission ratio of a 

-3 -4 

single interference filter is generally 10 to 10 

When this type of probe reflectance measurements were attempted, 

we discovered, unfortunately, that enough direct ruby light or fluorescence 

produced by the ruby light was received at the detector to overwhelm the 
o 

6328 A probe laser signal, A light-tight holder to accommodate a sequence of 
alternate filters and apertures was constructed to prevent any light from 
reaching the detector by paths through the edges or around the filters. Even 
when the number of filters in the stack was increased to nine to overcome 
possible bleaching in the first few filters, a direct pulse was still detected 
from the ruby laser. We concluded that mechanisms such as bleaching or 
fluorescence produced by the high-power ruby pulse prevented adequate 
discrimination by a stack of interference filters. 

Liquid solutions were also examined. The smallest ratio of trans- 
o o 

mission at 6943 A to transmission at 6328 A achieved was 0.05. This was 
produced by the Q-switch solution of cryptocyanine in acetonitrile (Eastman 
A10220) but was far from what was required to provide adequate discrimina¬ 
tion between the two lasers, 

31 

We were able to carry out specular reflectance measurements using 

a He-Ne probe laser when a Jarrell-Ash 0. 5 m Ebert spectrometer with a 
o 

single 6328 A interference filter in front of the entrance slit was used to 
provide wavelength discrimination. Even with this arrangement a small 
light pulse was detected when the He-Ne probe laser was turned off. While 
the specularly reflected beam is collimated and can be focused into the 
acceptance angle of a spectrometer, the output beam from the integrating 
sphere has a 2tt steradian solid angle and too great a light loss to utilize a 
spectrometer in carrying out total reflectance probe-laser reflectance 
measurements. 

Alternate approaches are: 1) use a dye laser with greater power out- 

o 

put and somewhat greater separation from the 6943 A ruby laser, or 2) 
choose a probe laser whose wavelength is longer than the ruby such as a 
cw YAG at 1.06|jl. Then cut-on transmission filters would provide much 
easier discrimination against the higher-power but now shorter wavelength 
ruby laser. One problem with this second approach is that the probe laser 
beam would not be visible and there would be much greater difficulty insuring 
that its focused spot was entirely within the focused ruby spot. 













VIII. CONCLUSIONS 


Thus far we may conclude that specular reflectance is a sensitive in¬ 
dicator of surface deformation. On the other hand, total reflectance measure¬ 
ments show that until the surface temperature of a metal target reaches the 
vicinity of the boiling point, no significant deviation is observed in the re¬ 
flectance from that given by a Drude free-electron model. These results are 

9 

in agreement with those of von Allmen et al. but contradict the conclusions 
of other workers that a decrease of - 50% in the total reflectance occurs at 
or possibly even before the melting point of a metal surface is reached. 

A most important result has been the preliminary verification of our 
initial postulate — that direct, real-time measurements of the target material '* 
optical properties with subnanosecond resolution can provide crucial, process- 
revealing signatures in following the interaction of an intense laser beam with 
a metal as the surface progressively undergoes heating, plastic deformation, 
slip, vaporization, ejection of liquid metal, plasma formation, etc. Measure¬ 
ment of the target's optical properties should be continued along with spectro¬ 
scopic measurement of the plasma generated. Electron and target specie 

32 

densities and temperatures can be compared with recent theoretical results. 

In this way our understanding of the laser-material interactions will improve 
and the high-power densities available from lasers can be more fully utilized 
in a wide variety of material processing applications. 


U 
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Stanley Bielaczy, Research Assistant, Microwave Research Institute. 

Robert DiFazio, Student Assistant, Department of Electrical Engineering 
and Microwave Research Institute. 






XII. LIST OF COUPLING ACTIVITIES WITH OTHER PEOPLE AND 
- GRWP T DOIN ’ GTEELATEP ' RESEARCH - 

A. Meetings Attended : 

1. International Conference on LASERS >7 8, Orlando, Florida, 
December 11 to 15, 1978. 

a. William T. Walter was chairman of Session I on Materials 
ana Interactions and presented a paper, "Reflectance Change 
of a Metal Surface During Intense Laser Irradiation, " with 
K. Park. 

b. Maurice Newstein presented a paper, "Kinetic Theory Treat¬ 
ment of Metal Evaporation Front, " with N. Solimene and 

J. Hammer. 

2. American Society for Metals Conference, "Applications of 
Lasers in Materials Processing, " Washington, D. C. , 

April 18-20, 1978. 

a. William T. Walter presented a paper, "Metal Reflectance 
Changes During Intense Laser Irradiation, " with K. Park. 

b. Maurice Newstein presented a paper, "Structure of the 
Evaporation Front From Boltzmann Equation, " with N. 
Solimene and J. Hammer. 

3. Society of Photo Optical Instrumentation Engineers, Annual 
Meeting, San Diego, California, August 27-30, 1979. 

a. William T. Walter presented an invited paper, "Reflectance 
Changes of Metals During Laser Irradiation, " during the 
seminar "Laser Application in Materials Processing." 

4. Metallurgical Society of the American Institute of Mining, 
Metallurgical and Petroleum Engineers, Fall Meeting, 
Milwaukee, Wisconsin, September 16-20, 1979. 

a. Kalinath Mukherjee presented a paper, "Nature of High 
Enegy Pulse Laser Induced Melting in Sn, Zn and Bi, " 
during the symposium, "Microstructure and Properties 
of Rapidly Solidified Non-Ferrous Metals." 

5. Eleventh International Quantum Electronics Conference, Boston, 
Massachusetts, June 23-26, 1980. 

a. Maurice Newstein presented a paper, "Laser Metal Inter¬ 
action in Vacuum. " 
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B. 


Seminar Speakers at the Polytechnic 


1. Dr, Alex Stein, Exxon Research Corporation, Linden, New Jersey, 
presented a seminar, "CF 4 Laser Research, " on April 7, 1978. 

2. Professor R. N. Sudan, Laboratory of Plasma Studies, Cornell 
University , Ithaca, New York, presented a seminar, "Two- 
Dimensional Turbulence in Low 0 Collisional Plasma, " on 

May 19, 197 8 . 

3. Dr. Richard L, Fork, Bell Telephone Laboratories, Holmdel, New 
Jersey, presented a seminar, "Dynamics of Photoexcited Gallium - 
Arsenide Band Edge Absorption with Subpicosecond Resolution, " on 
July 7, 1978. 

4. Dr. Sigrid McAfee, Bell Telephone Laboratories, Murray Hill, 

New Jersey, presented a seminar, "Nonlinear Absorption in 
Germanium on a Picosecond Time Scale, " on August 4, 1978. 

5. Dr. Raul Stem, Bell Telephone Laboratories, Murray Hill, New 
Jersey, presented a seminar, "New Results on Plasma Filamenta- 
tion, " on September 22, 1978. 

6 . Dr. N. J. Berg, Harry Diamond Laboratories, Adelphi, Maryland, 
presented a seminar, "Acousto-Optical Signal Processing, " on 
October 12, J 97 8 . 

7. Professor Shalom Raz, University of Houston-Texas, presented 
a seminar, "An Open Problem in Inverse Scattering Theory," 
on February 23, 1979. 

8 . Dr. M. Keskinen, Naval Research Laboratory, Washington, D.C., 
presented a seminar, "Numerical Simulation of the Rayleigh-Taylor 
Instability, " on March 9, 1979. 

9. Professcr C. Lee, Cornell University, Ithaca, New York, pre¬ 
sented a seminar, "Laser Annealing, * on March 30, i979. 

10. Dr. S, Basu, Air Force Geophysics Laboratory, presented a 
seminar, "Multi-Technique Observations of Equatorial Irregularities 
in the Night Time F-Region, " on May 4, 1979. 

11. Dr. E. McDonald, Naval Research Laboratory, presented a 
seminar, "Gradient-Drift Processes and Formation of Barium 
Cloud Striations, " on May 18, 1979. 

12. Dr. H. Eun, National Bureau of Standards, presented a seminar, 
"Gravity Waves and Atmospheric Species, " on May 25, 1979. 

13. Professor T. Shiosaki, Kyoto University, Japan, presented a 
seminar, "ZnO Thin Film for Bulk and SAW Applications, " on 
August 3, 1979. 

14. Professor I. Senitzky, The Technion, Israel, presented a seminar, 
"Cooperative Relaxation of a Coherently Driven Three Level System, " 
on October 5, 1979. 
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15. Dr. Hyatt Gibbs, Bell Laboratories, Murray Kill, New Jersey, 
presented a seminar, "Optical Bistability, ■ on January 11, 1980. 

16. Professor Joel DuBow, Colorado State University, presented a 
seminar, "Automatic Network Analysis of Materials and Devices, " 
on February 11, 1980. 

17. Dr. M. Lax, Bell Laboratories, Murray Hill, New Jersey, pre¬ 
sented a seminar, "Unstable Resonators in a Medium with Gain, * 
on February 15, 1980. 

18. Dr. S. McCall, Bell Laboratories, Murray Hill, New Jersey, 
presented a seminar, "Theory of Optical Bistability, " on 
February 22, 1980. 

19. Dr. Paul Rabinowitz, Exxon Research Laboratories, Linden, 

New Jersey, presented a seminar, "Controllable Pulse Com¬ 
pression in Raman Lasers, " on June 6, 1980. 
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